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METHOD AND APPARATUS FOR QUALITY OF SERVICE DETERMINATION 

This application is a continuation-in-part of U.S. application number 10/185,131, entitled 
"Catprobe", filed June 28, 2002. 

PRIORITY 

[0001] The present patent application claims priority to U.S provisional patent 

application serial nos. 60/457,902, entitled "Wirelesscat," filed on March 27, 2003 and 
60/533,186, entitled, "Catprobe," filed on December 29, 2003, and are hereby incorporated by 
reference herein. 

FIELD OF THE INVENTION 

[0002] The present invention relates generally to IP communication networks and more 

particularly to a quality of service estimation method for communication paths on IP networks. 

BACKGROUND OF THE INVENTION 

[0003] The "best-effort" nature of the Internet makes the QoS (Quality of Service) 

perceived by end users unpredictable and sometimes largely varying. Fast, accurate and efficient 
tools for estimating QoS performance of IP networks are gaining importance in the networking 
community. This is because such information can be used to maintain expected user and 
provider service under the varying conditions inherent to packet networks, especially the Internet. 
Specific applications include congestion control, real-time streaming and two-way 
communication, QoS verification, server selection and network administration. 
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[0004] QoS estimation can be broadly classified into two categories: passive monitoring 

and active monitoring. The passive monitoring approach has the advantage of not injecting 
additional probing traffic into the network. It observes the network as it is, meaning that the 
measurements are an assessment of true network behavior since this latter is not disturbed by 
probing traffic intended for those measurements. 

[0005] The monitoring can take different levels of granularity depending on the degree of 

processing, storage and resources available. Packet monitoring, for example, allows observation 
of packet-by-packet information such as packet delay variation, packet size distribution, and 
throughput between host pairs. A higher level, with less overhead, can be achieved by flow level 
measurements to record the total number of bytes transferred, the flow start and finish time, 
among others. 

[0006] The main advantage of passive probing techniques is that they do not introduce a 

load on the network they monitor, which also means they do not distort the network traffic and 
therefore produce realistic estimates. However, their handicap is that they rely on existing traffic, 
which is not guaranteed to have desired characteristics for certain measurements. Bottleneck 
bandwidth measurement techniques, for example, require a certain packet size distribution and 
inter-packet departure rate often not met. Traffic monitoring consists in passively observing 
traffic characteristics for the purpose of inferring network performance. 
[0007] SNMP (Simple Network Management Protocol) and RMON (Remote 

Monitoring) are the most widely adopted standards for passive monitoring and typically consist 
of management agents or probes installed at various network elements (hosts, routers, switches), 
a MIB (management information base) containing collected data from the agents, a management 
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station or console, which collects the information from the probes, and a protocol for the 
exchange of information between stations and probes. MIBs comprise several groups such as 
statistics collected at the physical or IP layers for packet sizes, CRC errors, and so forth. Traffic 
monitoring with administrative control requires the transfer of collected information from agents 
to consoles, thus placing a burden on the network being monitored. Sometimes, sampling of data 
in MIBs can be used to reduce the amount of traffic exchanged. 

[0008] Active monitoring obtains an inference of network QoS by sending probes across 

the network and observing the treatment they receive in terms of delay in delivery to the 
destination, variability in that delay and loss. A large variety of such tools exist to estimate 
performance in terms of delay, jitter, packet loss, and bandwidth. They generally either use 
Internet Control Message Protocol (ICMP) error messaging capabilities or packet dispersion 
techniques. 

[0009] Link capacity estimation in bits/sec has traditionally been achieved through the 

use of packet dispersion techniques, which consist of the successive transmission of groups of 
two or more packets. The concept is that packets from the same group will queue one after 
another at the bottleneck link of the path. With the absence of large interference from competing 
traffic (i.e. from other sources), the dispersion (i.e. the difference in packet arrival times at the 
receiver) will be inversely proportional to the bandwidth of the bottleneck. Examples of tools 
using this approach include Nettimer, Pathrate and Packet Bunch Mode (PBM). 
[0010] Another approach for capacity estimation builds on the ICMP Time Exceeded 

message option. Pathchar, the precursor of this technique, performs measurements by sending 
packets with increasing IP Time-to-Live (TTL) values, thereby forcing routers along the path to 
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send back ICMP error messages and revealing themselves. A measurement of round-trip delays 
to successive hops on the end-to-end path thereby leads to per hop delay estimation. Pathchar 
also adds the feature of varying packet sizes for each TTL value, thus inferring link capacity as 
the slope inverse of the line connecting minimum observed delay for each packet size. Other 
tools exist like Pchar and Clink but they build on the same concepts as Pathchar. 
[0011] Cprobe and Pipechar were tools proposed to estimate available bandwidth on a 

path. These tools use long packet train dispersion and assume that dispersion of such trains is 
inversely proportional to the rate available for transmission at the bottleneck hop, i.e., available 
bandwidth. Recent research, however, has shown that the inverse of such dispersions does not in 
fact measure available bandwidth but another parameter referred to as ADR (Asymptotic Data 
Rate). 

[0012] Another tool, Delphi, assumes Internet paths can be modeled by a single queue, 

which makes it perform badly in the presence of large queuing delays at several links on the path 
or when the bottleneck capacity and bottleneck available bandwidth links are located at different 
links. 

[0013] At the present time, there is only one known tool that is capable of estimating 

available bandwidth, which is Pathload. Pathload operates by performing several iterations, 
varying the transmission rate at each and observing delay variation at the receiver, with a goal of 
finding the maximum rate that does not cause delay increase. Pathload builds on the simple 
principle that the observed end-to-end delay at the receiver increases when the transmission rate 
at the source exceeds the available bandwidth on the path; a realistic observation knowing that 
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traffic injected onto the path faster than the bottleneck can service will cause queue build up at 
that hop, hence increasing queuing delay and delay altogether. 
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SUMMARY OF THE INVENTION 

[0014] A method and apparatus for estimating QoS in a network. In one embodiment, 

the method comprises probing an end-to-end path to identify addresses of all hops on the end-to- 
end path, generating at least one time-stamp request packet, transmitting the at least one time- 
stamp request packet to at least one hop on the end-to-end path, generating a time-stamp in 
response to the time-stamp request packet with the hop, and processing the time-stamp to 
produce at least one QoS estimate. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] The present invention will be understood more fully from the detailed description 

given below and from the accompanying drawings of various embodiments of the invention, 
which, however, should not be taken to limit the invention to the specific embodiments, but are 
for explanation and understanding only. 

[0016] Figure 1 illustrates an example of an IP network. 

[0017] Figure 2 is a block diagram of router packet processing. 

[0018] Figure 3 illustrates an exemplary varying of queue size with incoming rate. 

[0019] Figure 4 is a flow diagram of one embodiment of a process for generating QoS 

estimates. 

[0020] Figure 5 illustrates an exemplary ICMP time-stamp request processing at a 

desired pair of hops. 

[0021] Figures 6 and 7 show a model depicting the behavior of packets arriving at an 

access router (AR), traveling through a corresponding Access Point (AP), traversing a wireless 
link and reaching a remote terminal. 

[0022] Figure 8 illustrates the relationship between variables regarding access points to 

remote terminal communication. 

[0023] Figure 9 illustrates a Bernoulli model. 

[0024] Figure 10 illustrates a Gilbert model. 

[0025] Figure 11 is a graph depicting average and standard deviation of burst lengths. 

[0026] Figure 12 illustrates a three-state Markov model. 

[0027] Figure 13 illustrates a four-state Markov model. 
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DETAILED DESCRIPTION OF THE PRESENT INVENTION 

[0028] A method and apparatus for producing QoS estimates for communication paths in 

IP networks are disclosed. The communication paths may be between two given routers or hosts 
on an IP network. The QoS estimation technique may be used on a regular basis to poll a desired 
transmission path. One embodiment of a method of producing estimates comprises first 
identifying network addresses of all routers or hosts on the end-to-end path of interest and then 
transmitting pairs of time-stamp requests to each router or host on the end-to-end path. Finally, 
the time-stamp requests that are returned by the routers or hosts are processed to produce QoS 
estimates. For the purpose of the present invention, the term "hop" will be used collectively to 
refer to routers and hosts along the end-to-end path between respective terminals or nodes. 
[0029] In one embodiment, the technique to obtain addresses of hops on the path uses a 

Traceroute application. In addition, in one embodiment, time-stamp requests are obtained using 
ICMP time-stamp request/reply options. After using Traceroute to determine the addresses of all 
hops on the desired path, a group of ICMP time-stamp requests is sent from the source at one end 
of the path to each router on the desired path. In one embodiment, the group of ICMP time- 
stamp requests comprise five ICPM time-stamp requests. The first three packets of each group 
are used to determine initial conditions for the measurement to the corresponding host or router 
(i.e., hop on the path). The remaining two probing packets of the group of five are used to obtain 
an estimate of the total time spent by an ICMP request packet at the corresponding hop. 
[0030] Another pair of similar ICMP time-stamp request packets is sent to each hop to 

similarly estimate the total time spent by an ICMP request packet at each hop. These may be sent 
at regular intervals. Such successive measurements may be used to infer queuing delay variation 
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at each hop. In one embodiment, the estimates are formed at the originating source or node. In 
another embodiment, such estimates are formed at the receiver or destination node. 
[0031] In one embodiment, the total service time at each hop and the queuing delay 

variation are used to estimate the utilization in the direction of the path at each hop, which refers 
to the fraction of free transmission resources not used by competing traffic at the link. The 
estimates that are formed are also used to estimate the available bandwidth of each link in the 
direction of the path, which refers to the rate the link can further sustain without the occurrence 
of congestion and queue build up at the link. In addition, the estimates are used to estimate the 
rate of incoming traffic from competing sources to the link. 

[0032] In one embodiment, the technique described herein is performed by an application 

that resides at one node only that originates the QoS measurement. This node performs the 
processing required to generate QoS estimates based on collected time-stamp values from the 
hops. In another embodiment, this application resides at both source and receiver on the path of 
interest. In that case, the receiver may be in charge of processing the collected time-stamps and 
processing them. One difference in this embodiment is that ICMP packets sent to routers have an 
address spoofed to the destination such that replies are sent to the receiver instead of the source, 
which would have originated the ICMP requests. 

[0033] Further advantages of the present invention will be apparent from the following 

description, reference being made to the accompanying drawings wherein preferred embodiments 
of the invention are clearly illustrated. 

[0034] In the following description, numerous details are set forth to provide a more 

thorough explanation of the present invention. It will be apparent, however, to one skilled in the 
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art, that the present invention may be practiced without these specific details. In other instances, 
well-known structures and devices are shown in block diagram form, rather than in detail, in 
order to avoid obscuring the present invention. 

[0035] Some portions of the detailed descriptions which follow are presented in terms of 

algorithms and symbolic representations of operations on data bits within a computer memory. 
These algorithmic descriptions and representations are the means used by those skilled in the 
data processing arts to most effectively convey the substance of their work to others skilled in the 
art. An algorithm is here, and generally, conceived to be a self-consistent sequence of steps 
leading to a desired result. The steps are those requiring physical manipulations of physical 
quantities. Usually, though not necessarily, these quantities take the form of electrical or 
magnetic signals capable of being stored, transferred, combined, compared, and otherwise 
manipulated. It has proven convenient at times, principally for reasons of common usage, to 
refer to these signals as bits, values, elements, symbols, characters, terms, numbers, or the like. 
[0036] It should be borne in mind, however, that all of these and similar terms are to be 

associated with the appropriate physical quantities and are merely convenient labels applied to 
these quantities. Unless specifically stated otherwise as apparent from the following discussion, 
it is appreciated that throughout the description, discussions utilizing terms such as "processing" 
or "computing" or "calculating" or "determining" or "displaying" or the like, refer to the action 
and processes of a computer system, or similar electronic computing device, that manipulates and 
transforms data represented as physical (electronic) quantities within the computer system's 
registers and memories into other data similarly represented as physical quantities within the 
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computer system memories or registers or other such information storage, transmission or display 
devices. 

[0037] The present invention also relates to apparatus for performing the operations 

herein. This apparatus may be specially constructed for the required purposes, or it may 
comprise a general purpose computer selectively activated or reconfigured by a computer 
program stored in the computer. Such a computer program may be stored in a computer readable 
storage medium, such as, but is not limited to, any type of disk including floppy disks, optical 
disks, CD-ROMs, and magnetic-optical disks, read-only memories (ROMs), random access 
memories (RAMs), EPROMs, EEPROMs, magnetic or optical cards, or any type of media 
suitable for storing electronic instructions, and each coupled to a computer system bus. 
[0038] The algorithms and displays presented herein are not inherently related to any 

particular computer or other apparatus. Various general purpose systems may be used with 
programs in accordance with the teachings herein, or it may prove convenient to construct more 
specialized apparatus to perform the required method steps. The required structure for a variety 
of these systems will appear from the description below. In addition, the present invention is not 
described with reference to any particular programming language. It will be appreciated that a 
variety of programming languages may be used to implement the teachings of the invention as 
described herein. 

[0039] A machine-readable medium includes any mechanism for storing or transmitting 

information in a form readable by a machine (e.g., a computer). For example, a machine- 
readable medium includes read only memory ("ROM"); random access memory ("RAM"); 
magnetic disk storage media; optical storage media; flash memory devices; electrical, optical, 
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acoustical or other form of propagated signals (e.g., carrier waves, infrared signals, digital 
signals, etc.); etc. 



Overview 

[0040] Figure 1 is a block diagram of one embodiment of an IP network. Referring to 

Figure 1, an IP network 10 is generally illustrated that includes hosts or routers 12 and at least 
two terminals or nodes communication with each other through respective hosts or routers 12 on 
IP network 10. Information is communicated through the IP network 10 as packets or datagrams. 
The terms "packet" and "datagram" are used interchangeably herein to refer to a data package 
having a header. 

[0041] The collection of hosts or routers 12 and terminals or nodes at each end and the 

links connecting them constitute a path. When transmitting, a packet or datagram will make 
several hops in order to reach its destination. These hops may be across hosts or routers 12 and, 
as such, as used herein, the term "hop" should be construed to include both hosts and routers 12. 
The transmitting node or terminal can be referred to as a path source (or source, in short) and the 
receiving node or terminal can be referred to as a path destination (or destination, in short). As 
illustrated, source terminals can either be remote terminals 14 or fixed terminals 16. Likewise, 
destination terminals can be either remote terminals 18 or fixed terminals 20. 
[0042] In one embodiment, fixed terminals 16 and 20 directly connect through a fixed 

link to a respective router 12, remote terminals 14 and 18 connect through a wireless link 22 to 
an access point 24. An access point may include a wireless communication tower 26 and a router 
12. Alternatively, an access point 24 may include a wireless local area network (WLAN) 

006655.P029X 13 



compliant access point 26, such as an IEEE 802.1 1 access point with our without an associated 
router 12. 

[0043] Remote terminals 14 and 18 are illustrated as wireless telephones, but those 

skilled in the art would recognize that other remote terminals such as personal digital assistants 
(PDA), laptops, or the like might be used on the wireless IP network 10. As such, the depiction 
of a wireless telephone in Figure 1 should be viewed in an illustrative sense and not as a 
limitation of the present invention. As further illustrated in Figure 1, each router 12 is connected 
to at least one other router 12 to enable IP packet transmission among terminals across the IP 
network 10. 

WiredCat 

[0044] Figure 2 is a block diagram of one embodiment of IP packet processing in a router 

12. Referring to Figure 2, IP packets arrive at multiple ingress line cards 30 and enter an IP 
processing stage 32. Packets are then typically placed in an input queue 34 or an output queue 
36, depending on the specific implementation of the router 12. A switch fabric 38 directs packets 
to an appropriate egress line card 40. Queuing can be at either the input to switch fabric 38 or the 
output (although most modern routers implement input queuing). However, it always occurs 
after IP processing 32. As set forth in greater detail below, in the case of use of the ICMP time- 
stamp option, the time-stamp retrieved by the ICMP time-stamp request records the time before 
entry into input queues 34 or output queues 36. The IP processing time is small since packets are 
usually handled at wire speed at that stage so that waiting time and the actual processing time for 
IP processing stage 32 are negligible compared to the queuing in switch fabric 38 and the 
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transmission time at egress line cards 40. Queuing of a packet after IP processing 32 is due to the 
waiting time of other packets to reach their turn for transmission. Queuing latency directly 
relates to the size of the queue upon arrival of the packet under consideration and the 
transmission speed of the link; the packet needs to wait for the queue ahead to empty before it 
can be transmitted. 

[0045] In the discussion set forth below, the following notations are used and refer to the 

following items: t n or n in short is the continuous time at which an event is recorded; Qi(n) is the 
size in bits of the queue in link i at time n; Bi the capacity of link i in bits/sec; Ri(n) is the sum of 
all incoming flow rates to the queue in link i in bits/sec at time n; qi(n) is the queuing delay in 
router i at time n; Ui(n) is the utilization of link i at time n; and m is the packet size in bits. 
[0046] For a packet entering the queue at time n, the size of the queue ahead of it is Qi(n). 

The queuing delay q£n) of that packet is the time required to empty the queue in link i at rate Bi 
(i.e., the queuing delay of the packet), which can be represented as: 

qiin) = ~B~ (1) 

The rate Bi is referred to as capacity, the total bit forwarding speed of a link. In other words, it is 
the transmission rate of the link. This value is therefore fixed for a given router 12 and network 
wire configuration. 

[0047] Figure 3 depicts the change in queue size (hence, queuing time) with the variation 

of incoming flow rates. As long as /?, is less than the speed Bi at which the link can serialize bits 
on output linecard 40, the queue size will not increase. However, because of the bursty nature of 



006655.P029X 



15 



Internet traffic, it is possible that /?, may exceed the capacity 5„ thereby causing queue build up. 
For the purpose of the present invention, the utilization of link i at time n is defined by: 

u i (n) = R i (n)/B i (2) 

Referring to Figure 3, the following equation can be obtained: 

Q i (n) = max{ £ [J% (r) - B t ](t r - f M ), 0} (3) 

r=l 

or: 

Qi (n) = max{ [R t (n) - ](f n - ^ ) + ft (n - 1), 0} (4 ) 
where, for simplicity, it is assumed that upon network initialization at time zero, the queue at all 
links is empty, hence Qi(0)=0, Vi. 

[0048] Note that the time interval (t r _i, t r ) is not set to unity since this interval will be 

varied later. Using equations (1) and (3): 

max{i[^(r)-^](? r -r^), 0} 

* {n) = — ~ 5 < 5 > 



Using equations (1) and (4): 

, iW=mx{ !Mlh^Vd) t , ;( „. 1)i0| (6) 

Substituting equation (2) in (6) results in: 

q { (n) = max{ [u t (n) - \}{t n - t n _ x ) + q { (n - 1), 0} (7) 

Hence: 



006655.P029X 



16 



iij (n) = max {— + 1, 0} (8) 

[0049] The traversing time Tj(n) is the time needed for a chunk of m bits to traverse a 

link at time n, namely: 

Ti(n) = j- + qi (n) (9) 

The quatity m/Bi is also referred to as transmission time of the link; traversing time as defined 
herein is therefore the sum of transmission and queuing delays while neglecting propagation and 
processing delays. Using equation (7) in (9): 
m 

T i (») = —+ max { [ Ui (n) - 1] x (t n - t n _ x ) + q t (n - 1), 0} 
B i 

Hence: 



*,= (10) 

T; (n) - max { [u, (it) - \](t n - t n _ x ) + q. (n - 1), 0} 



Then, substituting equation (2) in (10): 



Ri T t (n) - max{ [«,(») - l](t n - t n _, ) + q 9 (n - 1),0} ° ° 

[0050] Available bandwidth for a new flow is the rate a link can support such that the 

queue size (i.e., queuing delay) does not increase. In other words, available bandwidth is the 
fraction of the capacity B of the link that is not used by competing flows with sum of rates R: 
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AvBw i (n) = B i -R i (n) (12) 
Substituting equation (2) in (12) yields: 
Avfiw l (n) = fi l [l-w i (n)] 

W/ (n)<l (13) 
AvBwrfn) is the available bandwidth of link i at time n. 

[0051] The available bandwidth metric can be used to assess the margin of additional 

input flows the path can take before congestion starts building up, hence leading to increased 
delay and jitter. Bandwidth requirements for voice coders, for example, are in terms of available 
bandwidth. 

[0052] Throughput is usually measured by observing the time required to get a certain 

amount of bits across a link or path. As such, throughput is the ratio of bits to total time to get 
them through. Throughput of link i as measured at the receiver can therefore be written as 
follows: 

k 

ThrBw t (w, k) = — ( 14) 

B i 

where ThrBwi(n } k) is throughput in bits/sec in link i at time n, k the size in bits of the data 
packets observed, which is also referred to as a data chunk (possibly composed of packets of 
different sizes), B\ the link capacity and q x the queuing delay of the k bits chunk at the link. 
[0053] Over a path comprised of several links, and using equation (14) the throughput 

equation becomes: 
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k 1 
ThrBw nnth (n, k) = n -* 

^ ie f ath ThrB Wi (n,k) 

The denominator on the right-hand side of equation (15) actually corresponds to the delay needed 
to traverse a link and path respectively. Throughput is therefore a good indicator of delay 
conditions in the network. As such, throughput does not in itself describe capacity or available 
bandwidth of the links or path as defined herein. It does relate to those metrics, however, in the 
following way: 

[0054] Equations (2) and (5) in the continuous domain also give: 

q L (n) = max { J [w,- (r ) - l]dr , 0} 
'0 

Then, using equation (13): 

qi (n) = max{- J -ilZrfr, 0} 

'0 B i 

Finally, this result in equation (14) gives: 

B: 

ThrBwi (/!,*) = max { l - , B t } ^ 

1--J AvBwi(r)dr 

[0055] Figure 4 is a general flow diagram of one embodiment of a process for generating 

QoS estimates. The process is performed by processing logic that may comprise hardware 
(circuitry, dedicated logic, etc), software (such as is run on a general purpose computer system or 
a dedicated machine), or a combination of both. The essence of the process consists of 
estimating the variation in queuing delay at successive hops. Then, equations (2), (10), (1 1) and 
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(12) above are used to deduce the capacity B, available bandwidth AvBw, sum of incoming 
competing flow rates R and utilization w. 

[0056] As depicted in Figure 4, the process begins by processing logic identifying 

network addresses of all routers 12 (hops) along the path from path source 14 and 16 to path 
destination 18 and 20 (processing block 40). In one embodiment, this is achieved using a 
Traceroute application. Next, processing logic obtains a first estimate of the queuing delay at 
each hop on the path (processing logic 42). Processing blocks 40 and 42 form an initialization or 
bootstrap phase 44 of the processor. From then on, pairs of ICMP time-stamp requests are sent 
to each router 12 on a regular interval. This allows monitoring queuing delay variation 46 at 
routers 12 using the equations described above (processing block 48). In one embodiment, each 
ICMP time-stamp request produces one estimate for each router 12 of the QoS parameters 
capacity, available bandwidth and throughput, which is illustrated as QoS estimates 50 in Figure 
4. In addition, link utilization and size of competing traffic is estimated. Processing blocks 46, 
48 and 50 form a measurement phase (52). 

[0057] In one embodiment, measurement begins by invoking Traceroute application to 

destination sources 18 and 20 to identify all the hops on the end-to-end path. Once the list of 
hops is obtained, pairs of ICMP time-stamp requests are used to obtain the variation delay at each 
of the hops. Below is a discussion of how this is achieved for one hop, the procedure is identical 
for the remaining hops on the path. 

[0058] Figure 5 depicts one embodiment of the processing of a pair of ICMP time-stamp 

requests 60 and 62 at consecutive routers 12. Upon arrival at a router 64, packet 60 receives a 
time-stamp from an IP processing unit 66. Then, in one embodiment, packet 60 returns to its 
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originator through a return queue 68 of router 64. Packet 62 continues through a forward queue 
70 through a link 72 and arrives at the IP processing unit 66 of router 74 where it receives a time- 
stamp. Then, in one embodiment, packet 62 returns to its originator through return queue 68 of 
the router 74. 

[0059] To obtain the queuing delay variation at the link from router 64 to the router 74 

depicted in Figure 5, two ICMP time-stamp requests are sent one after another with no time 
difference between the transmission of the last bit of the first packet and that of the first bit of the 
second packet. The ICMP request of the first packet is directed towards the first router 64, while 
that of the second packet is sent to the second router 74. Packets of the same pair are identified 
using an Identifier or Sequence Number field in the ICMP packet. 

[0060] The time-stamp at the first router 64 of the first packet of the pair is denoted by 

©i 1 ; and that of the second packet of the same pair at the second router 74 by 0 2 i+i- For the 
purpose of identifying different measurements on a time-scale and without loss of generality, it is 
assumed that the time-stamps reported by ICMP time-stamp requests of the same pair are both 
recorded at the time that the second packet of the pair receives its time-stamp. This is equivalent 
to the notation Q\ } Q) and 0 2 i+i(j)> meaning that both time-stamps are collected at the same time 
instant j. The following relation gives the difference between these two time-stamps: 

O) - ©/ (7) = f + 4,- O) + d^ M (j) 

where m is the size of the ICMP time-stamp request packet (namely 40 bytes; 20 bytes for each 
of ICMP and IP), q£j) is the queuing delay of the packet pair j waiting to be transmitted from the 
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first router 64 to the second router 74, and di^i(j) is the propagation delay between the two 
routers. 

[0061] In one embodiment, one assumption that is made is that the time-stamp of the first 

packet not only denotes the time at which the packet is leaving the IP processing stage at router 
64, but also that of the second packet at the same stage. This assumption requires the two ICMP 
packets to arrive after the other at router 64, i.e., no interfering traffic between the two packets 
until they reach router 64 to be queried. 

[0062] Neglecting propagation and processing delay, the equation obtained is: 

©? + iO>e!(;)=^-+^(j) (17) 

Equation (17) above is the same as equation (9), i.e., it defines the traversing time across the link 
72 from router 64 to router 74, taking into account the transmission and queuing times at that 
link and neglecting propagation and processing delays. 

[0063] After the transmission of the first pair, the sender waits for a user-defined delay a 

before sending the next pair. In turn, the measurements collected are 0}(y'+l) and of+iO' + l) . 
Therefore, the queuing delay variation is deduced as follows: 

[®M U + 1) " ©?+i O)] + [©) O) " ©/ U + 1)] = Qi U + 1) " Qi O) ( 1 8) 

Now using equation (8): 

u i (y + l) = m ax { t«CJ*P-«</)] + i,0) 
[QmU+V-GmU)] 
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where U|(j+1) is the utilization at link i 72 between router 64 and router 74 as observed between 
the two measurement pairs j and j+1 at time tj, respectively. In addition: 



of(; + i)-0fO) = a+X A ^^) 

[0064] From equation above, ^ (j) is the queuing delay variation happening at all 



\<k<i 



upstreaming nodes before measured node /. As a result, the queueing delay set 

{ Aq k (j), k e [1,/]} can only be solved iteratively. At first, the queuing delay at first node, 

Aq x (j) , is solved by 

A 9l (7) = 0f(j + l)-0 1 2 a)-^. 
[0065] Next, assume, for any k <= N node, Aq k (j) is known, Aq N+] (j) is calculated 

according to 

Now, the capacity can be determined using equation (17): 



m 



" i = ef(j)-e}(j)- qi u) (20) 

The unknown variable in equation (20) is the queuing delay at time tj in router 64. However, it is 
known that: 



006655.P029X 



9i 0) = 9, 0 - 1) + 0-1. 7) (2 1 ) 

where Aq,{j-l,j) is the variation in the queuing delay at router 64 as observed between packet 
pairs j-1 andj. Hence: 

?«0)= £ A^(r-l,r) + ^(0) (22 ) 

r=l 

[0066] In equation (22), 0,{0J is an estimate of the queuing delay encountered by packets 

of the first pair sent at the start of the measurement. An estimate of the first queuing delay qtfO) 
at router 64 is obtained at the very start of the measurement process; first, by sending three 
consecutive ICMP time-stamp requests, the first one to router 64 and the next two to router 74. 
This will lead to the following: 

0? + i(O)-0? +1 (O) = - m 



5,(0) 

Hence, an estimate of the capacity Bi at time 0 of the link from the first router 64 to the second 
router 74 is obtained: 

B i (0) = ~1 2 (23) 

e? +1 (O)-0? +1 (O) 

[0067] The estimate of the capacity Bi will later be refined as the measurements progress 

as described below. However, the equation (17) at time 0 can also be written: 

e? +1 (0)-©J(0)=-^-+^.(0) 

which, when plugged in to equation (23): 

q t (0) = 2 X 0? +1 (0) - ©j (0) - & M (0) (24) 
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This initial estimate of the queuing delay at link i can be used throughout the measurement 
process in equation (22). An estimate of the queuing delay can also be repeated at every 
measurement event through the same procedure of sending three consecutive ICMP packets to 
obtain a better estimate for use in equation (22). 

[0068] Once an estimate of q£j) is obtained, it is used in equation (20) to estimate link 

capacity Hence, B t can be estimated using equations (20), (22) and (24). Now, the sum of all 
incoming competing flows /?, can be estimated using equations (2), (19), (20), (22) and (24). 
This is achieved by first estimating capacity B iy then using equation (19) to estimate utilization 
and equation (2) to find /?,-. In addition, available bandwidth AvBwt can be estimated by using 
equations (2), (12), (19), (20), (22) and (24) through the same process. 

[0069] In one embodiment, all routers 12 are probed regularly and at the same frequency. 

In another embodiment, certain routers 12 can be probed more or less often based on queuing 
delay variation rate of those routers 12. 

[0070] In one embodiment, the path hop addresses are identified first before ICMP pair 

transmission begins. In another embodiment, such pair transmission begins for each hop as soon 
as that hop is identified. 

[0071] In one embodiment, ICMP time-stamp requests are the standard size as defined 

for ICMP protocol. In another embodiment, dummy data can be added in the IP payload field 
thereby increasing transmission time at the router, a useful feature for very fast links. 
[0072] In one embodiment, path identification and time-stamp requests are sent from 

either source or destination hosts. In another embodiment, different time-stamp requests and 
path identification can be sent from either source or receiver or other hosts on the network 
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simultaneously. In addition, processing of measurements can be achieved at any, some or all 
hosts on the network including source and destination hosts. 

[0073] A further extension to the technique described herein includes avoiding the 

transmission of time-stamp request pairs beyond the slowest link on a path. This is to avoid 
dispersion between time-stamps as they arrive at the link of interest. To achieve this, probing 
takes place from either source or destination such that the bottleneck link is avoided. In one 
embodiment, observing the rate of variation of the different links identifies the bottleneck link. 
Another extension to the method is to further refine the precision of measurements by taking into 
account propagation delay on the links. This is achieved using a propagation delay estimate. 

An Alternative Wireless Link Embodiment - WirelessCAT 

[0074] A method and apparatus for estimating an end-to-end path packet loss rate in a 

network is described. In one embodiment, the method includes probing an end-to-end path to 
identify queue capacity of all routers on the end-to-end path, transmitting a ping request packet to 
a hop on the path, and processing a ping reply to produce a QoS estimate. In another 
embodiment, a method and operations for calculating a QoS value for a wireless link of an end- 
to-end network path includes measuring an average packet loss value for a wired portion of the 
network, and calculating a packet loss value for the wireless link as a function of the packet loss 
value for the wired portion of the network. The discussion of the present invention set forth 
below describes an embodiment for use with a wireless link 22, and is referred to herein as 
WirelessCAT. 
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[0075] Referring to Figure 1 , methods that assess network congestion across IP network 

10 (such as by observing queuing delay variation at routers 12 along a desired path allow the 
estimation of available bandwidth AvBw. Although these techniques are used with wired 
networks, where congestion is solely due to queue build up at routers 12, in a wireless link 22, 
however, packet congestion may be due to several other factors, in addition to queuing 
congestion. These include contention for the medium in shared media technologies such as IEEE 
802.1 1, GSM, PCS, and the like. Furthermore, packet congestion in a wireless network may be 
due to retransmission delay caused by random packet losses on the wireless link. 

[0076] A method for measuring available wireless bandwidth is described for use with an 

illustrative wireless hop of an IEEE 802.1 lb wireless LAN. The method may be used with other 
technologies, and the use of an IEEE 802.1 lb wireless LAN (WLAN) for illustration herein is 
not intended to limit the invention claimed herein. The raw transmission rate at the physical 
layer in IEEE 802.1 lb changes dynamically between four rates: 1 Mbps, 2 Mbps, 5.5 Mbps and 

1 1 Mbps. The selected rate varies according to the perceived Frame Error Rate (FER). The 
thresholds on FER for switching to a higher or lower rate are implementation-specific. 
[0077] In one embodiment, the first step for estimating available wireless bandwidth is 
obtaining the current total raw capacity of the wireless link. After obtaining the current raw 
capacity of the wireless link, the next step is to estimate packet congestion to assess how much of 
the raw link capacity is available. 

[0078] The wireless link raw capacity BPHY at the physical layer may be determined by 

software generally provided to users by manufacturers of WLAN devices. For example, Cisco 
Systems provides Aironet "350" client utility software for Windows CE 2.1 1. The Aironet "350" 
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client utility software includes a site survey function whereby a client adapter may read the status 
of a wireless link every 250 msec. The Aironet "350" client utility software includes two 
operational modes: passive and active. In the passive mode, no additional RF traffic is generated 
by the utility, only actual traffic is monitored and a raw capacity BPHY is estimated. In the active 
mode, the software allows for setting the characteristics of active probing streams that will be 
transmitted to estimate performance. Such characteristics include the number of packets to be 
transmitted, packet size, packet inter-departure delay and raw transmission rate, which, in one 
embodiment, may be 1, 2, 5.5 or 1 1 Mbps. 

[0079] Packets traversing wireless links suffer from additional delays as compared to 

packets traversing wired links. On wired links, packets typically experience processing, 
transmission, congestion and propagation delays. On wireless links, packets may additionally 
experience contention and retransmission delays. As used herein, contention delay refers to the 
time required for a packet to "grab" a wireless channel, due to the shared nature of the medium. 
As used herein, retransmission delay refers to the time required to successfully transmit a packet 
across a lossy wireless link, which may require several attempts. Packet transmission on wireless 
links may be modeled as described below. 

[0080] Figures 6 and 7 show a model depicting the behavior of packets arriving at an 

access router (AR) 12, traveling through a corresponding Access Point (AP) 24, traversing the 
wireless link 22 and reaching the remote terminal 18. 

[0081] One of packets 126 destined for the remote terminal 18 first reaches the access 

router 12 where it experiences a queuing delay q CO ngestion in congestion queue 122 waiting for 
other of packets 126 ahead arriving from other routers 12 to be transmitted. Each of packets 126 
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emerge from the congestion queue 122 one after the other to enter into contention queue 124. 
When one of packets 126 comes to the head of the line in congestion queue 122, it experiences a 
contention delay waiting for the wireless link 22 to become available for transmission. The 
contention delay is caused by the shared nature of the transmission medium in technologies such 
as IEEE802.11b. 

[0082] Each of packets 126 becomes the head of the line in the contention queue 124, it 

is transmitted across wireless link 22. In case of failure, it reenters contention queue 124 from the 
tail. This is repeated until the packet is successfully received at remote terminal 18. However, 
there is a limit n on the number of times one of packet 126 is retransmitted beyond which the 
packet is simply dropped. In effect, a packet 126 entering the tail of congestion queue 122 on 
access router 12 experiences a total queuing delay q k until it reaches remote terminal 18. This 
delay is given by the equation following: 

0 ~ Q congestion + A t x +(A — 1)T+ X B\ ^^propagation (25) 



where: 



A k is the total number of transmissions for one packet k at the head of 



contention queue 124 at access router 12; 

Bf is the contention time for the i* retransmission of the packet k at 
the head of the contention queue 124 at access router 12; 

T is the timeout delay before a retransmission of the packet ahead is 
decided; this includes the roundtrip propagation delay and the transmission of the 
acknowledgement packet and processing at the remote terminal 18; 
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f * is the transmission delay of packet k over wireless link 22; 

k 

1 propagation 1S ^ propagation delay of packet k across wireless link 22; 

and 

Congestion ls the time s P ent the P acket ahead of packet k in both 
congestion queue 122 and contention queue 124 at access router 12. 

[0083] Equation (25) can be rewritten in term of averages and ignoring propagation delay 

(which is around O.S^isec for a typical 100m wireless link 22) as the following equation: 

t = ^congestion + A k t k +{A k -Y)T+A k B k (26) 

Since any packet will be retransmitted up to n times before it is dropped: 
A k <n, Vk (27) 
where n is typically equal to 7 in most implementations of DEEE802.1 lb. 

[0084] Two embodiments for estimating the number of transmissions A k are discussed 

herein. The first embodiment builds on information provided by the client utility of access point, 
while the second embodiment uses packet loss estimates. 

[0085] Turning to the first embodiment, existing utilities for remote terminals 18 or 

access points 24 typically accumulate information from Table 1 the moment the driver of the 
client utility of the access point is loaded or the moment wireless link 22 is established until the 
driver is unloaded or wireless link 22 is terminated. For example, the client utility for a Cisco 
Aironet "350" WLAN card accumulates the following data in real-time: 
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Table 1 



Access Point to Remote Terminal 



a=number of frames transmitted successfully. 



b=number of frames retransmitted. 



c=number of frames that exceed maximum number of retries. 



d=average contention time. 



The specific relationship among these variables is further illustrated in Figure 8. 

[0086] This accumulated data may be written into a text file. The text file may then be 

probed at specific time intervals T to calculate variations on each of the variables: Aa, Ab, Ac. 

The total number of packets handled during interval T is a+c. The total number of transmissions 

is: 



[a-(b-c)] \ +az + £A* <[a-(b-c)] + cn + (b-c)n = [a-(b-c)] + bn 

k=\ 

On average for one packet k: 

- k ^ bn + a -b + c (28) 
a -he 

and also: 

B k =d 

Using (26), the following is obtained: 

k k bn + a-b + c k ( bn + a-b + c (bn + a-b + c)d 

<1 ^4 congestion* ~ h H " 1K + " (29) 

a+c a+c a+c 



where: 



l*-Jt- (30) 
B PHY 
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where k is the size of an 802.1 lb frame. The timeout delay r is implementation specific, and 
^congestion iS obtained using WiredCAT congestion estimation described above. 

[0087] The second embodiment for estimating the number of transmissions A will now 

be described. Let p be the probability of packet loss on wireless link 22. The number of 
transmissions experienced by packet A: is a random variable A k where the probabilities for the 
occurrence of each value of A k is shown in Table 2 below: 



A 


# transmissions 


#retransmissions 


#losses 


#success 


probability 


1 




0 


0 




p u (i-p) 


2 




1 


1 




p'(i-p) 


3 




2 


2 




p z (i-p) 


j 




j-1 


j-1 




p>'(i-p) 


n-1 




n-2 


n-2 




P nz u-P) 


n 




n-1 


n-1 


0 


p ni (i-p)+p n 



Table 2 

[0088] On average, the number of transmissions required for successful transmission of 

packet k across wireless link 22 is given by: 



A* =1 Lprobability(i) = Zip 1 \\-p) + np n 



\e A k - l ~ P Tin* I nn"-nn" ■ - ( " + + P 

i.e., A — i^ip +np — np H 



n+1 



P /=1 
Again referring to (26): 



p<\-pY 



-k -k n np n+2 -(n + \)p n+} +p k . n np tn+2 -(n + \)p n + l + p t . „ np n+2 -(n + \)p n+l + p J 

P(\-P) 2 Pd-P) 2 p{\~p) 2 
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[0089] The fixed terminals 16 of the IP network 10 transmits to remote terminal 18, 

packets 126 traverse both wired and wireless hops. Therefore, end-to-end packet loss monitoring 
is not sufficient to infer wireless packet loss. In this scenario, a packet loss observed from the 
endpoints may be due to either a loss on the wired link because of congestion, or due to a loss on 
the wireless link because of congestion or random loss and does not correspond to the packet loss 
of wireless link 22. 

[0090] A number of implementations may be used to infer wireless packet loss at the 

network layer. Two such embodiments estimate the loss on the wired hop and then infer 
behavior on the wireless link from knowledge of end-to-end packet loss through RTCP and 
packet loss models described below. Both of the embodiments infer wired packet loss via 
functionality implemented at the path endpoints, in this case fixed terminals 16 and remote 
terminal 18. These embodiments are performed by processing logic that may comprise hardware 
(e.g., circuitry, dedicated logic, etc.), software (such as is run on a general purpose computer 
system or a dedicated machine), or a combination of both. 

[0091] In one embodiment, processing logic starts by inferring loss on the wired link 

alone for a system where software or hardware for QoS estimation is only deployed at the path 
edges. There are proposals to estimate QoS by deploying dedicated software on a router 12 of an 
access point 24. However, at least one embodiment of the present invention eliminates the 
requirement for dedicated software on a router 12 of an access point 24. 
[0092] In another embodiment, processing logic infers wired packet loss at the network 

layer using ping packets. The essence of this embodiment is to trigger a response from router 12 
of access point 24 and infer wired hop packet loss without actually deploying specialized 
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software at the router 12. For that purpose, ping packets are interleaved within a voice packet 
stream. These pings target router 12 which is revealed by a Traceroute operation on the path. The 
rate of ping packets interleaving is preferentially set in such a way as to capture packet loss 
behavior while not exceeding a given overhead to the traffic. 

[0093] For example, sending one ping every twenty voice packets effectively results in an 

overhead of five percent in packet count. Since packet loss in wired networks occurs in bursts, 
observing the loss of ping packets may provide a good approximation of the overall voice 
packets loss rate. The time interval between ping values is set or varied in order to capture the 
burstiness of packet loss; as such, it is not a regular interval. The interval may vary in order to 
capture the behavior of a wired network packet loss model. Exemplary packet loss models are 
described below under the heading "Packet Loss Models." 

[0094] The second embodiment infers wired packet loss via additional functionality at 

the path endpoints. This additional functionality may be implemented as a third stage added to a 
wired links embodiments described above, in addition to the initialization and measurement stage 
described above. This new stage is dubbed "router buffer estimation." In one embodiment, the 
new stage operates as follows. With knowledge of available bandwidth on wired links, the 
sender transmits a stream of ICMP time-stamp pairs such as is described above at a rate higher 
than the estimated available bandwidth to each router 12 on the wired hop from fixed terminal 16 
to access point 24. The queuing delay at each router is inferred and monitored until ICMP 
packets start getting lost at the targeted router 12 due to queue overflow. For each router 12 
number /, the last observed queuing delay value q^i is recorded. Then, the queue capacity at 
each router 12 number i is estimated as being: 
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Qi =9max4 B i (31) 
where B t values are wired link i raw capacity that has been measured. This router buffer 

estimation stage may occur only once at the start of the measurement process. Thereafter, 

WiredCAT may keep repeating the measurement stage until the real-time measurement is ended. 

The path may be probed every T seconds. An estimate of wired packet loss pd during the probe 

interval T is derived, with the assumption of constant conditions during interval T t as follows: 

If Ri(T)<Bi then: 

Pd-0 during T. 

If Ri(T)>Bi then: 



R + r-B 
R + r 

Pd = 0 during Fi(T) 



Pd ~~ n during T-Fi(T); and 

K + Y 



So: 



where: 



Pd ~ j X R + r during T. (32) 



where Q is in bits, q is in seconds, B and R are in bits/second; R is the total background traffic 
entering the queue and r is the traffic rate of the application. FtfT) is the time needed to fill up the 
queue assuming R remains constant during T. 
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Packet Loss Models 

[0095] A method to combine wired and wireless packet loss pd and p w respectively to 

produce path packet loss p t will now be described. This expression provides a method for 
inferring p w after calculating pa and given a knowledge of p t using the methods discussed above. 
Good models for packet loss in IP networks include the Bernoulli random loss model, the 2-state 
Markov chain model (or Gilbert model) and the more complex k th order Markov chains model. 
Embodiment of the present invention makes use of the former two models, which are described 
in the following paragraphs. 

[0096] The Bernoulli model is a two-state Markov chain model that can be used to model 

random and independent packet losses. Figure 9 below illustrates a Bernoulli model where "0" 
represents a state of successful packet transmission and "1" represents a loss. The transition 
probability from state "0" to "1" is denoted by /?, and the transition probability from state "1" to 
"0" is denoted by q. The relation between q and p is as follows: 

= (34) 

[0097] The Gilbert model, on the other hand, is a two-state Markov chain model that may 

represent the burstiness of packet loss due to congestion in IP networks. Figure 10 illustrates a 
Gilbert model where "0" represents a state of successful packet transmission and "1" indicates a 
packet loss. The transition probability from state "0" to "1" is denoted by p y and the transition 
probability from state "1" to "0" is q. To incorporate the effect of burstiness, q and p do not sum 
to one as in the case of the Bernoulli model. Instead, the following equation holds: 
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where r is the total average packet loss. 

[0098] The average and standard deviation of burst lengths for the Gilbert model are 

illustrated in Figure 11. The burstiness of packet losses Bl relates to the parameters of the model 
by: 

Bl= £ nq{\-q) n - X (36) 

[0099] Next, path packet loss will be described in order to derive a relation between 

measured packet losses on the wired and wireless hops and the end-to-end loss that characterizes 
the entire path. The derivation of the model is divided into two steps. First packet loss models 
are derived separately for the wireless and wired hops. Then, these two models are combined 
into a single model. 

[00100] For the wired hops of the path, the Gilbert packet loss model is used. It is 
assumed that losses on the wireless hop are due to both congestion at access point 24 and to 
random losses on wireless link 22. Again, the Gilbert model is used to handle congestion at 
access point 24, while the Bernoulli model is used to account for random packet losses on 
wireless link 22. An expression is then derived to relate end-to-end packet loss p t on a path 
between endpoints 1 and 2 to the wired hop loss information. 

[00101] Packet losses on the wireless hop can be due to either congestion at the access 
point 24 or bit errors at wireless link 22. These two stages are modeled separately as a and b. 
Figure 12 illustrates a 3-state Markov model that represents the losses on the wireless hop for a 
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single packet experiencing both congestion at access point 24 and random loss on wireless link 
22. 

[00102] In Figure 12, p QtS is the conditional packet loss due to congestion, pb, s is the 
unconditional packet loss due to random effects, and q QtS reflects the unconditional packet loss 
due to congestion. Parameters relating to congestion are noted with the subscript notation a and 
parameters pertaining to random losses have the subscript notation b. All of the parameters have 
the subscript s to show their relation to wireless link 22. 

[00103] The representation in Figure 12 uses three states. A value of "0" indicates no loss 
and a value of "1" signifies a loss. State (a,b) reflects the conditions of transmission (loss or 
success) at the congestion stage or wireless link 22 stage, respectively. The state (1,1) is 
effectively the same as (1,0) since a packet lost on the wired hop will not continue to the wireless 
hop. The three-state model shown in Figure 12 may be completely described by the following 
transition probability matrix: 



Ps = 



(! " Pa,s )(! " Pb,s ) C 1 - Pa,s ) Pb lS Pa,s 
(! " Pa,s X 1 " Pb 9 s > ( J " Pa,s > Pb § s Pa,s 
Va,s (X " Pb,s ) Va,s C 1 - Pb,s ) <X-4a 9 s) 



(37) 



[00104] To obtain P 0 ,o , Poj , and P u , the system is solved as follows: 

Po.ol 1 





"p 

r 0,0 


1 




P o,. 




{ 






. p u 


J 


I 


P H-P 

r 0,0 T r 0, 



0,1 



V.J 



(38) 
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[00105] The end-to-end path average packet loss p t wireless link average packet loss is 
then given by: 

r s = \\ + P U (39) 
[00106] Equation 39 assumes that packet loss on separate hops is independent. The model 

does, however, take into account correlation in packet losses on the same link from packet to 

packet. 

[00107] The packet loss model for the entire path, including the wired and wireless hops, 
is obtained by combining the three state model of Figure 12 with the Gilbert model, as shown in 
Figure 13. In Figure 13, pd is the conditional packet loss on the wired hop and qj reflects the 
unconditional packet loss on the wired hop. Parameters with the subscript d relate to the wired 
hop and those with s to the wireless hop. Additionally: 

CCx=P 2 =iX-p d ){\-p a ,s)^-Pb,s) 

a 2 =0 l =(l-p d )(l-p as )p btS 
a 3 =fo=(l-p d )p atS 
a 4 =J3 4 =S 2 = p d 
S l =a-P d )(l-P a , s ) 

S 3=( l ~Pd)<la,sPb,s 

^=i}-Pd)q a ^-P b ,s) < 40 > 

<Pl=<ldQ--<la,s) 
9i = Wa,sPb,s 
<P4=<ld<laA X -Pb,s) 



[00108] The model shown in Figure 13 uses four states, where each state describes the 
conditions on the wired and wireless hops. A "0" value means no loss and a "1" signifies a loss. 
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State (d,a,b) reflects the conditions of a transmission as a loss or a success at the wired hop, 
congestion stage of the wireless hop, or wireless hop stage, respectively. The state (0,1,1) is 
effectively the same as (0,1,0) since a packet lost on the wired hop will not continue to the 
wireless hop. Similarly, states (1,1,1), (1,0,0), (1,1,0) and (1,0,1) are identical. 
[00109] The four-state model shown in Figure 13 may be completely described by the 
following transition probability matrix: 



Pt = 



a x a 2 a 3 a 4 

Pi Pi Pi Pa 

8 X S 2 S 3 S 4 

<P\ <Pi <h <P4 



(41) 



To obtain P 0 ,o,o , Po,o,u Poj,i , and P U j , the matrix is solved yielding: 

l 



^0,0,0 


1 


"p 

^0,0,0 


1 


p 


= Pt 


p 




^0,1,1 




p 




P.,u 


J 


_P,x. 


J 




+ P 0.U 





(42) 



Where the end-to-end path average packet loss p t is given by: 

Pt =P 0,0,1 + ^0,1,1+ P \X\ (43) 
[00110] Equation (43) assumes that packet loss on separate hops is independent. The 

model does, however, take into account correlation in packet losses on the same link from packet 

to packet. 
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[00111] 



The four-state Markov chain model shown in Figure 13 may be used to calculate 



p t from measured average packet losses on the wired and wireless hops rj and r s , respectively. To 
calculate p t using the model, the first step is to calculate pd, q& Pa,s* Pb,s and q QtS . A value for qt 
and q a , s may be obtained by maintaining a measurement of average packet burst length and using 
equation (36) and the graph shown in Figure 1 1, as described above. A starting value for both qj 
and tf^may be 0.9, corresponding to an average loss of 1.1 1 packets/burst and standard deviation 
of 0.35 packets/burst. (See Figure 1 1). These values have been shown to be the most probable 
for Internet, currently. As the Internet evolves, other starting values may be preferred. 
[00112] From equation (35), pd may be obtained as follows: 



[00113] The value of p a , s may be obtained by solving the following system of three 
equations and two unknowns { p a , s \ Pb.s } • 




(45) 



(1 " Pa,s )(! " Pb,s ) (1 " Pa,s )Pb,s 
Ps= (X- P a ,s)(X- Pb,s) ( l -Pa,s)Pb,s 



Pa,s 
Pa,s 



4a,sd- Pb,s) <laA l -Pb,s) 
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where r s is the measured average packet loss on the wireless link, q a>s is calculated as described 
above. Once the values of p& q& Pa,s, Pb.s and q a , s are known, equation (42) may be solved for 
Po.o.0, Po.o.h Po,u> and P UtL Finally, plugging the values for P 0 ,o,o> Po.oj , Pojj , and P U j 
plugging into equation (43) will yield the end-to-end path average packet loss p u 
[00114] Whereas many alterations and modifications of the present invention will no 
doubt become apparent to a person of ordinary skill in the art after having read the foregoing 
description, it is to be understood that any particular embodiment shown and described by way of 
illustration is in no way intended to be considered limiting. Therefore, references to details of 
various embodiments are not intended to limit the scope of the claims which in themselves recite 
only those features regarded as essential to the invention. 
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